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Abstract 

Boson-boson scattering and Higgs production in boson-boson fusion hold the 
key to electroweak symmetry breaking. In order to analyze these essential features 
of the Standard Model we have performed a partonic level study of all processes 
I1I2 —> Q'.j Qa at the LHC using the exact matrix elements at 0 (a® m ) provided 
by PHASE, a new MC generator. These processes include also three boson production 
and the purely electroweak contribution to tt production as well as all irreducible 
backgrounds. Kinematical cuts have been studied in order to enhance the VV 
scattering signal over background. PHASE has been compared with different Monte 
Carlo’s showing that a complete calculation is necessary for a correct description of 
the process. 
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1 Introduction 


The nature of Electro-Weak Symmetry Breaking (EWSB) will be investigated at the 
LHC 1 . The Standard Model (SM) provides the simplest and most economical explanation 
of the phenomenon and all its predictions have been verified with unprecedented accuracy. 
The only missing ingredient is the Higgs boson which is essential to the renormalizability 
of the theory. The fit of precision EW data to the SM currently gives an upper limit 
on the Higgs mass of about 200 GeV, depending on the preferred value of the top mass 
[5]. A heavier Higgs mass can however be accomodated by a number of models which 
in general would be distinguishable from the SM by improved precision data or by other 
observations at future colliders [6] 

In the SM the Higgs is also crucial to ensure that perturbative unitarity bounds are not 
violated in high energy reactions. Scattering processes between longitudinally polarized 
vector bosons are particularly sensitive in this regard and without a Higgs they violate 
unitarity at about one TeV. 

It should be noted that the Goldstone theorem and the Higgs mechanism do not 
require the existence of elementary scalars. It is conceiveable and widely discussed in the 
literature that bound states are responsible for EWSB. 

Since unitarity is essentially a statement of conservation of total probability it cannot 
be violated in Nature. Violation of perturbative unitarity implies that the SM becomes 
a strongly interacting theory at high energy. The nature of the interaction between 
longitudinally polarized vector bosons and the Higgs mass, or possibly the absence of 
the Higgs particle, are strongly related: if a relatively light Higgs exists than the Vl’s 
are weakly coupled, while they are strongly interacting if the Higgs mass is large or the 
Higgs is nonexistent [7]. In the latter case, by analogy with low energy QCD, which 
can be expressed by exactly the same formalism which describes the Higgs sector in the 
SM, or adopting one of the many schemes for turning perturbative scattering amplitudes 
into amplitudes which satisfy by construction the unitarity constraints, one is led to 
expect the presence of resonances in WlWl scattering. Unfortunately the mass, spin 
and even number of these resonances are not uniquely determined [8]. If, on the other 
hand, a Higgs particle is discovered it will nonetheless be necessary to verify that indeed 
longitudinally polarized vector bosons are weakly coupled at high energy by studying 
boson boson scattering in full detail. If the Higgs mass is not far from the present 
experimental lower limits it will take several years of data taking before enough statistics 
could be accumulated for a reliable discovery. In the meantime high energy vector boson 
scattering could provide an alternative method for investigating the existence of a light 
Higgs scalar, providing a strong incentive to probe harder the low mass range if no sign 
of strong scattering is detected. 

At the LHC no beam of on shell EW bosons will be available. Incoming quarks 
will emit spacelike virtual bosons which will then scatter among themselves and finally 
decay. These processes have been scrutinized since a long time [9, 10]. All previous 
studies of boson boson scattering at high energy hadron colliders have resorted to some 
approximation, either the Equivalent Vector Boson Approximation (EVBA) [11], or a 
production times decay approach, supplementing a calculation of 

qiq-2 -> q-mViV-2 (l) 

1 Detailed reviews and extensive bibliographies can be found in Refs.[l, 2, 3, 4] 
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processes with the on shell decay of the two vector bosons. This is not surprising since 
the number of reactions required by a full six fermion calculation is very large and often 
approximate results are adequate for a first discussion. There are however issues that 
cannot be tackled without a full six fermion calculation like exact spin correlations between 
the decays of different heavy particles, the effect of the non resonant background, the 
relevance of off mass shellness of boson decays, the question of interferences between 
different subamplitudes. Without a complete calculation it will be impossible to determine 
the accuracy of approximate results. 

Recently a full fledged six fermion Monte Carlo has become available [12], based on 
the methods of Refs. [13, 14], which in its present version describes at 0(a® m ), using 
exact matrix elements, all processes of the form 

PP -> q^ 2 -> (2) 

(where q % stands for a generic (anti)quark) which can take place at the LHC. This provides 
a complete description of VV scattering and its EW irreducible background. Eq.(2) 
includes about one thousand different reactions. A complete classification of all processes 
can be found in [12], Preliminary results have been presented in Refs. [15]. 

Since in addition to VlVl scattering many other subprocesses are in general present 
in the full set of diagrams, it is not a trivial task to separate boson boson scattering 
from the EW irreducible background. In practice one has to deal also with other types of 
background to which QCD interactions contribute, but which however do not include any 
boson boson scattering term. We will refer to these processes as QCD background even 
though in general they will be a mixture of QCD and EW interactions, tt production 
and decay is particularly relevant and dangerous. In this paper we are neglecting QCD 
backgrounds, whose calculation is in progress, and concentrate on 0{oP em ) processes. It is 
clear that obtaining a good signal over EW irreducible background ratio is a prerequisite 
to any attempt at dealing with the QCD one. 

In the absence of firm predictions in the strong scattering regime, trying to gauge 
the possibilities of discovering signals of new physics at the LHC requires the somewhat 
arbitrary definition of a model of VlVl scattering beyond the boundaries of the SM. 
Some of these models predict the formation of spectacular resonances which will be easily 
detected. For some other set of parameters in the models only rather small effects are 
expected [8]. 

The simplest approach is to consider the SM in the presence of a very heavy Higgs. 
While this entails the violation of perturbative unitarity, the linear rise of the cross sec¬ 
tion with s, the invariant mass squared in the hard scattering, will be swamped by the 
decrease of the parton luminosities at large momentum fractions and, as a consequence, 
will be particularly challenging to detect. At the LHC, for Mh >10 TeV, all Born di¬ 
agrams with Higgs propagators become completely negligible in the LInitary gauge, and 
the expectations for all processes in Eq.(2) reduce to those in the Mh —> oo limit. In this 
paper we will compare this minimalistic definition of physics beyond the Standard Model 
with the predictions of the SM for Higgs masses within the reach of the LHC. 
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Figure 1: Vector boson fusion processes 



Figure 2: Non fusion and non doubly resonant two vector boson production. 




Figure 3: Higgs boson production via vector boson fusion and Higgstrahlung. 



Figure 4: Electroweak tt and single top production. 
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Figure 5: Three vector boson production. 


E(lepton)> 20 GeV 
Pt (lepton) > 10 GeV 
|//(lepton) | < 3 
E(quarks)> 20 GeV 
Pt (quarks) > 10 GeV 
1 77 (quark) | < 6.5 
M(qq)> 20 GeV 

Table 1: Standard acceptance cuts applied in all results. Any pair of colored fermions 
must have mass larger than 20 GeV. Here lepton refers to / ± only. 


2 The six fermion final state processes 

Boson-boson scattering and Higgs production in boson-boson fusion produce intermediate 
states with two bosons and two quarks as shown in Fig. 1. In this study we have only 
considered final states in which one W boson decays leptonically and the other (either Z 
or W) hadronically, which is regarded as the best channel for probing these processes at 
the LHC [3]. If both bosons decay hadronically the signal cannot be distinguished from 
the QCD non resonant background whose cross section is much larger. Final states where 
both vectors decay leptonically have a smaller rate and have been left for future studies. 

Once vector bosons are decayed we have six fermion final states. Since the V’s emitted 
by the initial state particles are spacelikc, we are forced by gauge invariance to include all 
diagrams in which final state V’s are emitted directly by the fermion lines as in the left 
part of Fig. 2. When the finite width of the EW boson is properly taken into account and 
the outgoing vector bosons are allowed to be off mass shell, it becomes necessary, again 
because of gauge invariance, to consider the full set of diagrams described by Eq.( 2 ). As 
a consequence, in addition to the diagrams which are related to the process we would like 
to measure, VV fusion, shown in Fig. 1, there will be diagrams in which a pair of V’s 
are produced without undergoing VV scattering, as presented in Fig. 2. Furthermore, 
diagrams related to Higgs production via Higgstrahlung will also be present as shown in 
Fig. 3, as well as diagrams which can be interpreted as tt EW production or as single top 
production as shown in Fig. 4. Finally diagrams describing three vector boson production 
which include Triple Gauge Coupling and Quartic Gauge Coupling will contribute as well, 
as shown in Fig. 5, since they produce the same kind of six fermion final states. 

Depending on the flavour of the quarks in Eq.(2) the various subprocesses will con¬ 
tribute and interfere to a different degree. All processes will be experimentally indistin- 
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guishable, apart from heavy quark tagging, and will have to be studied simultaneously. 

A number of samples of events representative of all possible processes in Eq.(2) have 
been produced with PHASE. In order to comply with the acceptance of the CMS detector 
and with the CMS trigger requirements, the cuts in Tab. 1 have been applied. We have 
used the CTEQ5L [16] PDF set with scale 

Q~ = M\v + g y^, (3) 

where pxi denotes the transverse momentum of the i-th final state particle. 

2.1 Physical sub-processes 

As already mentioned many subprocesses (i.e. WW ->WW, ZW —»ZW, ZZ ->WW, tt) 
will in general contribute to a specific six fermion reaction. It is impossible to separate 
and compute individually the cross section due to a single subprocess, since there are 
large interference effects between the different contributions. 

We can however select all complete 2^6 processes which include a specific set of 
subdiagrams. For instance, ZW —>ZW with on shell bosons is described by 4 Feynman 
diagrams. These same diagrams, with all external vector bosons connected to a fermion 
line, constitute the ZW —>ZW set of 2^6 diagrams. Several sets can contribute to a 
single process and therefore the same process can appear in different groups. The upper 
part of Fig. 6 shows the invariant mass distribution of the two most central quarks (when 
ordered in pseudorapidity rj) ,the lepton and the neutrino for the reactions which contain 
the different subprocesses as well as the distribution for the complete set of processes. We 
assumed M(H)=500 GeV. For this figure, in addition to the standard cuts decribed in 
Tab. 1, the following cuts have been applied. The pseudorapidity of the two most central 
quarks should satisfy | rj (central-quark)) < 3, the pseudorapidity of the two most forward- 
backward quarks are constrained to \rj (forward-quark)| > 1 and the transverse mass of 
the lepton neutrino pair should be smaller than 100 GeV. These three last requirements 
have been released for the studies presented in the following sections. 

It should be clear that the total cross section in Fig. 6 is smaller than the sum of the 
cross sections for the various groups. Notice that the Higgs peak is present in the ZW 
—>ZW curve. This is due to processes that in addition to the ZW —>ZW set of diagrams 
include also diagrams describing Higgs production, e.g. uu —> uuudp~V. 

The groups comprising single top and tt diagrams have a large cross section. An in¬ 
variant mass analysis reveals that the processes in these two groups are indeed dominated 
by top production. The lower part of Fig. 6 shows the same distributions after subtraction 
of single top and tt and with the additional requirement that the two most central quarks 
have invariant mass between 60 and 100 GeV. 

Top candidates are identified requiring a b-quark and two other quarks in the final 
state of the right flavour combination to be produced in a W decay, with a total invariant 
mass between 160 and 190 GeV. Analogously, events in which the muon, the neutrino and 
a b quark have an invariant mass between 160 and 190 GeV are rejected. These vetoes 
reduce drastically the EW top background and produce a much sharper Higgs peak. 
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Invariant Mass of 2 central q, u and v 



Invariant Mass of 2 central q, u and v 


all processes 
W W -> W W and C.C. 



Figure 6: Invariant mass distribution of the lepton, neutrino and the two most central 
quarks, for different sets of processes. In the upper part with the set of cuts described in 
the text, in the lower part after vetoing top and tt production and requiring the vector 
bosons to be close to their mass shell. 
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Invariant Mass of bosons 



Figure 7: The invariant mass of the two most central quark in pseudorapidity (red/dashed) 
and for the quarks that gives the best nominal vector boson mass (black/dash-dotted) 
with PHASE. The invariant mass of the lepton and neutrino pair is shown in blue/continous 
for comparison. 


2.2 Identification of the hadronically decaying vector boson 

In the LHC environment it is not obvious how to identify the jets produced by the hadronic 
decay of an EW vector boson among the four which are present at generator level in the 
processes we are interested in. Fig. 7 shows the invariant mass of the candidate vector 
boson selected with two different procedures. In the first one final state quarks are ordered 
in pseudorapidity. The quarks with the largest and smallest rj are identified as tag quarks, 
while the remaining two, which will be referred to as “most central” quarks, are identified 
as the quarks from the W or Z decay 2 . In the second procedure the two quarks among 
the four that have the mass closest to a vector boson mass are identified as the quarks 
from the decay while the remaining two are identified as tag quarks. As a comparison the 
invariant mass of the lepton and neutrino pair is also shown. It is clear from the figure 
that the second choice reproduces much better the lepton-neutrino mass distribution, 
giving higher efficiency and purity at generator level, and will be adopted in the rest of 
this paper. 

2 Notice that these two quarks are actually those with smaller \r)\ only in the center of mass of the 
hard reaction and not in the laboratory frame. However the ordering in pseudorapidity is preserved by 
boosts along the beam line and the selection of the two jets with intermediate rapidity will be the same 
in the two reference frames 













3 The VV-fusion signal 

In order to isolate the VV fusion process from all other six fermion final states and 
investigate EWSB, different kinematical cuts have been applied to the simulated events. 

The main background we want to subtract at parton level is single top and tt produc¬ 
tion. These events have been identified as in Sub-Sect. 2.1 and rejected. 

The invariant mass of the muon and the neutrino has to reconstruct the mass of a 
W, and is required to be in the range M w ± 10 GeV. In VV fusion an additional W 
or a Z decaying hadronically is expected to be present. Therefore events are required to 
contain two quarks with the correct flavours to be produced in W or Z decay, with an 
invariant mass of ± 10 GeV around the central value of the appropriate EW bosons. If 
more than one combination of two quarks satisfies these requirements, the one closest to 
the corresponding central mass value is selected. This combination will in the following 
be assumed to originate from the decay of an EW vector boson. 

In a second step, in order to reject events which can be identified with the production 
of three vector bosons, the flavour content and the invariant mass of the two remaining 
quarks is compared with a W and a Z. If compatible within 10 GeV with either, the 
event is rejected. This happens in about 2% of the cases. The events satisfying all these 
constraints will constitute the “signal” sample. 

These requirements are not fully realistic: no flavour information will be available for 
light quarks, b’s will be tagged only in the central part of the detector and the invariant 
mass of the Iv system will not be directly available. Our aim however is to define a 
“signal”, in the same spirit as CC03 was adopted as such at LEP 2, that is a pseudovariable 
which could be used to compare the results from the different collaborations. The signal 
is not necessarily directly observable but it should be possible to relate it via Monte Carlo 
to measurable quantities. If such a definition is to be useful it must correspond as closely 
as possible to the process which needs to be studied and the Monte Carlo corrections 
must be small. At this stage we want to isolate as much as possible the VV fusion signal 
from all other production channels, while keeping reasonably close to the experimental 
practice and taking into account the full set of diagrams required by gauge invariance. It 
becomes then possible to analyze the differences between the VV fusion signal sample and 
its intrinsic background. This also provides some preliminary experience at the generator 
level which could guide more realistic and complete studies. 

In Tab. 2 the cross sections obtained for each of the described cuts are shown. The 
corresponding VV mass distributions are displayed in Fig. 8 and 9. 

Fig. 8 shows how the invariant mass distributions for two reconstructed vector bosons 
is modified when the cuts described above are sequentially imposed, for the case where 
the Higgs boson exists and its mass is 500 GeV and for the no-Higgs case. The no-Higgs 
amplitude is the amplitude in the limit Mh —■► oo. In the Unitary gauge, for the processes 
under investigation, this limit amounts to neglecting all diagrams with Higgs propagators. 
The black line refers to the full sample, the red one shows the effects of antitagging on 
the top. The blue line corresponds to imposing all the mentioned cuts and antitagging 
on the presence of three vector bosons. Fig. 9 compares the reconstructed vector boson 
invariant mass distribution after all cuts for the two Higgs cases 3 . 

Since our definition of signal begins with the identification of the jets which reconstruct 

3 The black and red line of Fig. 9 do not reproduce the results of Fig. 6 because of a different set of 
cuts 
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Invariant Mass of 2 central quarks, n and v 


Invariant Mass of 2 central quarks, n and v 




Figure 8: Distributions of the invariant mass of the two central quarks, the muon and the 
neutrino for M(H)=500 GeV (left) and for the no-Higgs case (right). The black line refers 
to the full sample. The red one shows the effects of antitagging on the top. The blue line 
corresponds to imposing all the mentioned cuts and antitagging on the presence of three 
vector bosons. 


| Invariant Mass of VW | 



Figure 9: Distributions of the invariant mass of the two reconstructed vector bosons for 
M(H)=500 GeV (red) and for the no-Higgs case (blue) after all cuts.In this plot the vector 
boson that decays hadronically is reconstructed using the two quarks that give its best 
mass. 
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| Eta of quarks - quarks from boson decay 



Energy of quarks 



E [GeV] 


Delta Eta between quarks 




p [GeV] 


Figure 10: Distributions of the pseudorapidity 77 , the difference Arj, the energy and trans¬ 
verse momentum of the quarks. In red (full line) for the quarks from the decay of the 
vector boson and in blue (dotted line) for the tag quarks. All plots refer to M(H)=500 
GeV. The distributions for the no-Higgs case lead to the same conclusions. 


the EW boson rather than with the identification of the tag jets, we compare in some 
detail the kinematics of the quarks from a vector boson decay with the kinematics of the 
spectator ones. Fig. 10 shows that a good separation in phase space between the two 
pairs of jet is maintained. The distribution of the pseudorapidity 77 of the quarks coming 
from the decay of the vector boson and of the tag quarks shows that the quarks from 
the V decay are, as expected, rather central (i.e. low 77 ) while the spectator quarks tend 
to go forward/backward. Consequently the difference in pseudorapidity A 77 between the 
quarks from the bosons will be smaller than the A 77 between the tag quarks. The quarks 
from the vector boson have less energy and smaller transverse momentum with respect 
to the tag quarks. 

It is interesting to look at possible differences in the kinematics of the VV-fusion signal 
with respect to the irreducible background. In Fig. 11 the “no-Higgs” case is chosen as an 
example, but there are no significant differences with the case of a massive Higgs. Only 
some variables, which are connected to the mass of the Higgs boson, show the presence 
of the resonance. In that case the figure for the finite mass Higgs is shown. The total 
invariant mass of the six fermions in the final state is presented for the signal and for 
the full sample: the signal tends to have a very large final six fermion mass. The muon 
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from the signal has a larger pr than the one from the background, and the same applies 
to the spectator quarks. In the background the energy of the tag quarks are peaked at 
low values, although there is a long tail at high energy as for the signal. The transverse 
momentum of the W that decays leptonically is shown for the signal events only and for 
the total. We can clearly see the presence of the Higgs resonance around p t ~ Mjj/ 2. 
The difference of the rf s of the tag quarks is also shown: the signal tends to have a larger 
At] with respect to the background. 



no Higgs 

500 GeV 

170 GeV 

u (pb) 

perc. 

u (pb) 

perc. 

u (pb) 

perc. 

total 

0.689 

100 % 

0.718 

100 % 

1.003 

100 % 

signal 

0.158 

23% 

0.184 

26% 

0.409 

41% 

top 

0.495 

72% 

0.494 

69% 

0.496 

49% 

non VV resonant 

0.020 

3% 

0.023 

3% 

0.040 

4% 

three bosons 

0.016 

2 % 

0.017 

2 % 

0.057 

6 % 


Table 2: Cross sections and percentages of events for different Higgs masses 


4 PHASE versus other Monte Carlo’s 

In this section we compare PHASE with PYTHIA [17] and MADEVENT [18]. Since PHASE is 
the only dedicated code that uses an exact calculation for the 2 ^6 processes, it has 
been used as a reference. The PYTHIA Monte Carlo employs the Equivalent Vector Boson 
Approximation, and includes only longitudinally polarized vector bosons (processes 71- 
77). The cross section in PYTHIA depends strongly on the cut applied to the px of the 
outgoing W’s in their center of mass. This cut, stored in CKIN(3), is unavoidable since 
it eliminates the divergence due to the photon exchange diagram in the on shell WW 
scattering. Notice that this divergence is an artifact of the on shell projection of the 
2^2 amplitude. In the exact six fermion matrix element the photon divergence is 
absent. For the present comparison CKIN(3) was set at 50 GeV. 

Fig. 12 shows the comparison between PHASE and PYTHIA for the Mh = 500 GeV and 
for the case where no Higgs boson is produced (in PYTHIA this can be simulated by asking 
for a very high Higgs mass). The same signal selection cuts, including those in Tab. 1 
have been applied to both PHASE and PYTHIA events. For a Higgs of 500 GeV one finds 
a reasonable agreement around the resonance but a large discrepancy in the rest of the 
mass spectrum. For very large Higgs masses the two results are in sharp disagreement 
over most of the spectrum, but for the largest VV invariant masses. 

As a check, in Fig. 13 the contributions to the cross section due to the different 
polarizations in ud —> udcspV are shown for the two different scenarios. Since it is 
impossible to separate the polarizations of the vector bosons in the full 2-^6 matrix 
element, we have used a dedicated code that keeps only the diagrams in which two W 
bosons are produced. This corresponds to convoluting the matrix element for ud —> 
udW + W~ with the matrix elements for the polarized decay of the W’s and with their 
Breit-Wigner distributions. This approach violates gauge invariance, but since we are 
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I Six fermions mass | 




| P T of the W decaying leptonically | 



Muon p T 




Delta Eta between tag quarksH 



Figure 11: Differential cross section as a function of the invariant mass of the six fermions 
in the final state, the transverse momentum of the muon, the transverse momentum of the 
two tag quarks, the energy of the two forward quarks, the transverse momentum of the W 
decaying leptonically and the difference in pseudorapidity between the two tag quarks for 
all the events (black) and for the signal events (red). All plots refer to the no-Higgs case 
with the exception of the pt of the W which decays leptonically which is for M(H)=500 
GeV. 
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restricting the W’s to be close to their mass shell these effects are expected to be small. 
The cuts shown in Tab. 3 have been applied in this case, in order to enhance WW —> WW 
compared to ZZ —> WW. 

The TT cross section is essentially independent of the Higgs mass and it is the largest 
one, outside the peak region, for M H = 500 GeV. In the no-Higgs case LL and TT 
production contribute equally to the spectrum. An attempt to experimentally separate 
LL and TT polarizations has been done and is described in Sect. 6. 

For Mh = 500 GeV case the disagreement at high invariant mass is due to the fact 
that only LL is considered in PYTHIA, while the cross section for TT polarized vector 
bosons dominates. For the no-Higgs case, the cross section for LL and TT polarizations 
are of the same order. The agreement between PYTHIA and PHASE at high masses is 
probably accidental, the missing polarizations being compensated by the growth of the 
LL cross section in PYTHIA, which as the total center of mass energy increases gets larger 
contribution from the photon exchange diagram at small scattering angles. We have 
not pursued this issue any further since Fig. 12 makes it clear that PYTHIA is not an 
appropriate tool to investigate VV scattering. 


E(quark,lepton) > 20 GeV 
Pt (quark,lepton) > 10 GeV 
1 < r)(d ) < 5.5 
—1 > r)(u ) > —5.5 
70 GeV < M(cs,/x V) < 90 GeV 

Table 3: Acceptance cuts for ud —> udW + W~ —> udcsfiv in Fig. 13. The initial state 
u(d)-quark is in the z(—z) direction. 


MADEVENT in principle could simulate all processes using exact matrix elements. In 
practice it is far too slow to do so and is used in the “production times decay” approx¬ 
imation, i.e. computing only the processes in which two vector bosons are produced on 
their mass shell and then decaying them. Even if this procedure partially accounts for spin 
correlations, it neglects all correlations within decays. This approach overlooks a number 
of contributions and cannot be applied in the region below the VV threshold which will 
be actively investigated at the LHC in the quest for an intermediate mass range Higgs. 
Moreover, it leads to double counting in some regions of phase space. For instance the 
reaction 

ud —> udcsn~v (4) 

would be approximated by the sum of two processes 

ud —> udW + W~ —> udc§n~V ud —> W + csW~ — > udcsfi~v. (5) 

While this is a reasonable approximation over most of phase space since the ud pair from 
the first reaction will mainly be produced at large rapidity separation and the ud pair 
from the decay of the W + in the second reaction will be at small separation, and as a 
consequence interference effects will be small, it leads to an overestimate of the cross 
section in the region of three vector boson production where all final state particles are 
rather central. 
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Invariant Mass of VW 




Figure 12: Cross section as a function of the invariant mass of the two vector bosons for 
PHASE (red) and for Pythia (blue). The upper plot refers to a Higgs mass of 500 GeV, 
the lower one to the no-Higgs case. 
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We have compared the predictions of PHASE and those of MADEVENT for the process 
uu —> uuqq'nis with qq' = du + sc and for the sum of all processes with no b quarks 
(b’s require a separate generation in MADEVENT). In this case we have adopted the scale 
Q 2 = Mf since our standard choice Eq.(3) cannot be used in the production times decay 
approach. The signal has been selected as described in Sect. 3 therefore eliminating the 
region of three vector boson production which is not correctly described by MADEVENT. 
A number of distributions are shown in Fig. 14 and Fig. 15 respectively. Fig. 14 shows 
significant differences between the two results. The transverse momentum of the spectator 
quarks is harder in PHASE, as is the pseudorapidity of the reconstructed vector bosons, 
while the separation in pseudorapidity between the spectator quarks tends to be smaller. 
In PHASE the muon is more likely to fly in the direction of the W momentum than in 
MADEVENT. While none of these differences is dramatic their impact should be carefully 
assessed. 

These differences are also present in Fig. 15, only slightly diluted in the more inclusive 
sample. Fig. 15 presents also the distribution of the invariant mass of the two candidate 
vector bosons and of the Ar/ between the two most central quarks. At large invariant 
masses the two predictions coincide while close to the VV threshold MADEVENT predictions 
are about 7% larger. The difference in pseudorapidity between the two most central quarks 
is significantly harder in MADEVENT with differences of about 15% at large separation. 

It should be noted that the result of the comparison between the two codes is sensitive 
to the mass window around the nominal mass of the EW bosons which is included. In real 
simulations this will depend on the necessity of rejecting backgrounds and on the uncer¬ 
tainty with which invariant masses of two quarks and two leptons will be reconstructed. 
MADEVENT predictions show no dependence on this cut, since all W’s are produced exactly 
on mass shell. PHASE, on the contrary, takes into account the Breit-Wigner shape of the 
propagator squared and includes a large number of additional diagrams which become 
important when the vector boson are allowed to be off shell. 

5 Higgs production in PHASE 

Higgs production in VV fusion followed by Higgs decay to WW or ZZ is the second most 
abundant production channel over almost the full range of Higgs masses which will be 
explored at LHC. It is regarded as the channel with the highest statistical significance 
for an intermediate mass Higgs [19]. PHASE is capable of simulating Higgs production in 
VV fusion together with all its EW irreducible background for any Higgs mass and is 
particularly useful in the intermediate mass range where only one of the vector bosons 
can be approximately treated in a production times decay approach. For an intermediate 
mass Higgs the dilepton final state H —> WW (*) —> lulu is slightly favored with respect 
to the H —> WW^ —> lujj channel because of the W + nj background which affects 
the latter. In both cases the main background comes from tt production followed in 
importance by EW WWjj production which is estimated to be about 10%, much larger 
than QCD WW jj production. The latter can be reduced with a central jet veto which 
does not affect the former.. The cross section after selection cuts is of the order of 14-4 fb 
and with such a small number of events it is important to have a precise simulation of the 
complete final state with full spin correlations which are complicated by the fact that in 
this range of Higgs masses one of the W is far off shell. 
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In Fig. 16 we present a number of kinematical distributions for M(H)=140 GeV with 
different selection procedures. The green line refers to identifying the two most central 
quarks in pseudorapidity as the quarks from W decay while the black histogram refers to 
selecting the two quarks among the four that have the mass closest to the W boson mass. 
The blue line is obtained selecting the quark pairs which have the correct flavour content 
to be produced in a W decay and which have the mass closest to the W mass. The first 
procedure, which is also the most robust from an experimental point of view, agrees well 
with the third one, which is more appropriate from a theoretical perspective but which 
can only be applied at generator level, over the full range of all variables. It should be 
pointed out that the actual ratio of signal over background depends cruciallly on the 
experimental resolution and on the set of additional selection cuts [19]. For reference, the 
invariant mass distribution in Fig. 16 is plotted with a bin size of 10 GeV. 


6 The high mass region 

An interesting physics possibility is to investigate whether there exist or not an elementary 
Higgs boson by measuring the VW cross section at large M(VW). Preliminary studies 
performed at CMS with a fast detector simulation [20] showed that a resolution of about 
10-15% on M(VV) up to 2 TeV is achievable with about 100 fb _1 . The rise of the cross 
section related to unitarity violation in the no-Higgs case is difficult to detect at the 
LHC, since the center-of-mass energy is still rather low and the decrease of the proton 
distribution functions at large x has the dominating effect. As we discussed in the previous 
section, if WlWl final states could be isolated, the difference in the cross section at high 
WW masses would be sizeable, since the WlWl cross section decreases much more rapidly 
in the presence of a Higgs particle than in the no-Higgs situation. In order to distinguish 
WlWl from WtWt we must exploit the different behaviour of the final state in the two 
cases. 

To this purpose kinematical distributions for M(H)=170 GeV, M(H)=500 GeV and 
the no-Higgs case have been compared for M(VV)>800 GeV since the cross section at 
large M(VV) for M(H)=170 GeV and M(H)=500 GeV is essentially due to transversely 
polarized vector bosons, while the cross section for the no-Higgs case is due to a mixture 
of the two polarizations as shown in Fig. 13. Fig. 17 shows that the distributions are quite 
insensitive to the value of the Higgs mass provided it is much smaller than the invariant 
mass of the VV system. This raises the interesting possibility of defining Standard Model 
predictions for high invariant mass production of VV pairs. These predictions will obvi¬ 
ously suffer from the usual PDF and scale uncertainties, which could however in principle 
be controlled by comparing with the cross section of some appropriate “standard candle” 
process. 

We have tried several sets of cuts and we believe that using Neural Network is the most 
effective way of increasing the separation between the no-Higgs case and the presence of 
a relatively light Higgs. Two samples of events with a high invariant mass VW pair, for 
M(H)=500 GeV and the no-Higgs case respectively, have been employed to train a Neural 
Network. All events satisfy the cuts in Tab. 1. A set of variables which discriminate 
between the two Higgs hypotheses have been used in the training. We have chosen four 
weakly correlated variables: the difference in pseudorapidity between the two bosons and 
between the two tag quarks, the transverse momentum of the tag quarks and the cosine 
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of the angle between the lepton and the W boson in the W center of mass system. These 
(and other) kinematical variables are shown in Fig. 17 for the no-Higgs case and for 
M(H)=500 GeV and M(H)=170 GeV. 



no-Higgs case 

m(Higgs) = 500 GeV 


a 

£=100/fe _1 

<7 

£=100/6~ 1 

ratio 

all events 

13.6 fb 

1360 ± 37 

11.6 fb 

1160 ± 34 

1.2 

NN >0.52 

3.93 fb 

393 ± 20 

2.70 fb 

270 ± 16 

1.5 

NN >0.54 

3.17 fb 

317 ± 18 

1.95 fb 

195 ± 14 

1.6 

NN >0.56 

2.67 fb 

267 ± 16 

1.47 fb 

147 ± 12 

1.8 

NN >0.58 

2.28 fb 

228 ± 15 

1.13 fb 

113 ± 11 

2.0 


Table 4: Integrated cross section for M(VW)>800 GeV and number of expected events 
after one year at high luminosity. 


The neural network, a multilayer perceptron with BFGS training, available in the 
ROOT package [21], takes these four variables as input, has two intermediate layers, with 
eight and four neurons respectively, and one output variable. The differential cross section 
da/dNN where NN is the neural network output variable is shown on the left hand side 
of Fig. 18 for two event samples corresponding to the different Higgs masses. Applying a 
cut on NN we can enhance the separation between the heavy and light Higgs case. The 
corresponding purity and efficiency is shown on the right hand side. 

In Tab. 4 the integrated cross section and the number of events for the two Higgs 
cases for M(VW)> 800 GeV are shown for different values of the cut. The corresponding 
ratio is also presented. The predictions for M(H)=500 GeV have a statistical uncertainty 
of about 10% for NN > 0.56, while the very heavy Higgs case, with all the necessary 
caveats, predicts a cross section which is larger by about a factor of two. This is quite 
encouraging and suggests that further study should be worthwhile in order to sharpen 
the SM predictions. 

In Fig. 19 we show the ratio 


f°° AM da no Higg S 

r 

iM cut aivi vw 


( 6 ) 


for different cuts on the neural network output variable as a function of M cut , the lower 
limit of the integration range on Myw- The difference between the two cross sections 
is concentrated at large values of NN and increases with increasing values of the output 
variable. With a fixed lower cut on NN the ratio between integrated cross sections in¬ 
creases with M cut . With M cut > 1 TeV and NN > 0.56 the cross section is about 2 fb for 
a very heavy Higgs, roughly twice the expected yield for a light Higgs scalar. On the right 
hand side of Fig. 19 we show the ratio 6 for different group of processes. The set which 
includes W ± W ± is the one with the largest separation, while the sets including ZZ and 
ZW scattering are less sensitive to the cut on NN. 
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7 Summary and future 

In this paper we have studied all q±q 2 —> q-^-uq^qd^ processes at order 0(a® m ) at the LHC 
using for the first time a full fledged six fermion Monte Carlo event generator. We have 
examined how simple kinematical cuts can be applied at generator level to extract the 
VV signal from the irreducible background. In the high mass region we have compared 
the case of a relatively light Higgs with the no-Higgs case as a guide to separate the LL 
component of VV scattering. Employing a neural network approach it seems possible to 
obtain a good separation of the two cases. 

Work is in progress to extend PHASE to cover qqZZ final states. The calculation of 
the full set of processes 2 —> 6 at O(ola^ m ) is under way. A phenomenological study of 
VV processes with both vectors decaying leptonically is in project. 
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ud -> ud W + W~ -> ud p v c s~ 



ud -> ud W + W~ -> ud 4 v c s~ 



M(WW) [GeV] 


Figure 13: Cross section as a function of the invariant mass of the two vector bosons for 
their different polarizations. The upper plot refers to a Higgs mass of 500 GeV, the lower 
one to the no-Higgs case. 
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Figure 14: Cross section for the processes uu —> uuqq'fiu with qq' = du + sc as a func¬ 
tion of the separation in pseudorapidity between the spectator quarks, their transverse 
momentum, the pseudorapidity of the reconstructed vector bosons and the angle of the 
muon with respect to the W direction of flight in the W rest frame for PHASE (black) and 
for MADEVENT (red) for a Higgs mass of 500 GeV. 
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Invariant Mass of 2 q with M(V), jj and v 


xIO ' 3 




Eta of the bosons | 



Delta Eta between the tag quarks 





Figure 15: Cross section for the full sample as a function of the invariant mass of the two 
candidate vector bosons, the A// between the two most central quarks, the Ar] between 
spectator quarks, their transverse momentum, the pseudorapidity of the reconstructed 
vector bosons and the angle of the muon with respect to the W direction of flight in the 
W rest frame for PHASE (black) and for MADEVENT (red) for a Higgs mass of 500 GeV. 
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W=2q with M(W) and flavour 






Figure 16: Distribution of the invariant mass Myy of the two candidate vector bosons, 
their energy, pseudorapidity and transverse momentum for a Higgs mass of 140 GeV. In 
the last three plots only events with 120 < Myy < 160 GeV are included. 
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| Invariant Mass of 2 q with M(V), u and v 





| Lepton angular distribution (W center of mass RF) | 



| Eta of the central quarks 



| Eta of the tag quarks - MH170 signal 



| P T of the tag quarks 

xIO ' 3 




Figure 17: The invariant mass of the two vector bosons, the pseudorapidity r] of the two 
central jets, the Ar/ of the two central jets, the r) of the forward quarks the A rj of the 
two forward quarks, the transverse momentum of the forward quarks, the Ar] of the two 
vector bosons, the cosine of the angle between the lepton and the W boson in the W boson 
rest frame for M(VW) > 800 GeV. In green (full line) for the no-Higgs case, in black 
(dashed) for M(H)=500 GeV and in red for M(H)=170 GeV. The last four variables are 
used as input to the neural net. 
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| Neural Network output 


| purity and efficiency VS NN output cut~| 




Figure 18: da/dNN where NN is the neural network output variable, in red (full line) 
for the no-Higgs case and in black (dashed) for M(H)=500 GeV on the LHS. The corre¬ 
sponding purity and efficiency on the RHS 



M( vw ) cu( [GeV] 


No Higgs / M(H) = 500 GeV for different subprocesses (NN Qut >0.58) 



Figure 19: On the left, the ratio of cross sections, integrated between M cut and infinity, 
as a function of M cut , for different values of the neural network variable. On the right,the 
ratio Eq.(6) for different groups of processes under the condition NN > 0.58 
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